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CORRELATION OF IRRADIATION DATA USING ACTIVATION FLUENCES 
AND IRRADIATION TEMPERATURE 
by John H. Lynch 
Lewis Research Center 

SUMMARY 

Nil-ductility transition temperature data for A302-B ferritic steel were analyzed by 
using multiple regression. Four independent variables were used. These were speci- 
men temperature during irradiation and time -integrated specific activations (activation 
fluences) of fast-, intermediate-, and thermal-neutron detectors. The results of these 
analyses were predictive equations having the increase in transition temperature as the 
dependent variable. For a reference point, a similar analysis was also performed 
using the flux integral above 1 MeV and the irradiatior temperature as independent var- 
iables. Also, the effects of excluding irradiation temperature from these analyses was 
studied. 

The results of these analyses showed that excellent correlation was achieved by 
using activation fluences and irradiation temperature as the independent variables; that 
is, 93 percent of the data variability was explained as contrasted with 36 percent when 
using the flux above 1 MeV. Irradiation temperature was found to account for from 
29 to 46 percent of the data variability depending on the model used. 

The activation fluence model was also compared with a selected damage function 
model. This comparison indicated that the activation fluence method was about as ac- 
curate (for the data studied) as the damage function method, and had several additional 
advantages over this method. 


INTRODUCTION 

Radiation damage and the development of correlations for predicting damage have 
been the subject of many investigations (refs. 1 to 6 are typical). The correlations 
generally are attempts at relating a change in the physical property of a material to 
time-integrated neutron flux above some selected energy, usually 1 MeV but frequently 



as low as 0. 01 MeV. Those correlations which assume a threshold energy for damage 
are generally recognized to be inadequate for many situations because of the assump- 
tions made in deriving and using them. The energy distribution of neutrons is usually 
assumed to be that of a fission spectrum; only neutrons above 1 MeV are assumed to 
contribute (and those equally) to the damage production; and the observed damage is as- 
sumed to depend only on the total exposure and not on the exposure rate (ref. 7). The 
first two assumptions are clearly not generally true, especially for fast reactor applica- 
tions. The third has some theoretical (ref. 6) and experimental (refs. 8 and 9) verifi- 
cation. 

A more fundamental approach is typified by that of McElroy, Dahl, and Serpan 
(ref. 10), in which the damage is treated as a cross section (damage function) which 
when multiplied by the neutron flux yields the damage in both integral and differential 
form. Damage function methods usually involve complex calculations (ref. 10). They 
require knowledge of the neutron flux in some form (differential or integral above some 
selected energy) both for correlating data and for using these correlations to predict 
damage. This introduces error in predicted values because of uncertainties in cross 
sections and in the spectrum determination and error propagation in the damage function 
generation. 

The work discussed in this report correlates damage with detector activations. In 
the discussion that follows, this method is shown to have several advantages for corre- 
lating damage when compared with selected integral flux and damage function methods. 


METHOD 


The activation fluence method is essentially a regression analysis (ref. 11) of data 
to obtain a polynomial in several variables, their second powers, and their interactions. 
The variables account for spectrum shape (fast, intermediate, and thermal) and the 
temperature at which the damage is produced. 

In general, the polynomial is 


Y - A 0 + AjXj + A 2 X 2 + . . . + A n X n + A 12 X x X 2 + AjgXjXg + . . . + A^X^ 

+ a 23 x 2 x 3 + a 24 x 2 x 4 + . . . + A n _ ljn x n _ lXn+ a u x 2 x + a 22 x 2 


+ • • • + A nn X n + e 


(1) 


The dependent variable Y is the predicted damage. The coefficients A are determined 
by using multiple linear regression and are best least- squares estimates assuming only 
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that the experimental error e is normally distributed with a finite variance and a mean 
of zero. The X values are functions of the activation fluences for selected neutron 
detectors (Xj, Xg, . . X n l ) and temperature X n . 

Using multiple regression as a means to implement correlation gives an unambig- 
uous and immediately useful model along with quantitative indications of the accuracy 
that can be expected when using this model. Also, statistical interpretations are pos- 
sible which may aid in determining which variables are important so that more basic in- 
vestigations can proceed on an efficient experimental basis. 

The activation fluence for any detector is defined by 

/*T S' °° 

a = y Q y Q cr(E)<p(E,t)dE dt (2) 

where T is the total irradiation time and the other symbols have their usual meanings. 
In practice, a set of detectors is selected that will cover the complete range of neutron 
energy. When feasible, these detectors are irradiated with the irradiation specimen; 
and their activations are used to obtain the activation fluences. Detectors having long 
half-lives are preferable because they reflect the irradiation history more accurately. 


TESTING THE METHOD 

In order to test this correlation approach using actual data, a literature search was 
made for a consistent set of data. A data set was selected (refs. 12 and 13) that was 
consistently defined and that had accompanying calculated neutron spectra. These data 
represented a wide range of spectra and were taken in several different test reactors. 
The data and tabulated spectra along with several helpful comments were provided by 
C. Z. Serpan of the Naval Research Laboratory, Washington, D.C. 

The selected parameter of interest was the change in nil-ductility transition tem- 
perature ANDT for A302-B ferritic steel. Of the data examined, this parameter had 
by far the most data points. 

Detectors were selected in an attempt to cover the complete range of neutron ener- 
gies with minimal response overlap. Consideration was also given to availability, ease 
of counting, performance at high neutron flux levels, and expected accuracy of various 
detector materials. 

1 59 60 

The detector set first selected was Co (n, y)Co , responding from ~0 to 6 eV; 

Cu °(n, y)Cu° , responding from ~0. 4 eV to 0. 2 MeV when cadmium -covered; 

*In a fission spectrum 90 percent of the activation that would occur in these detec- 
tors is contained within the energy limits shown. 
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Np^^(n, f)N^, responding from ~0. 5 MeV to 5 MeV; and Ti^®(n, p)Sc^®, responding 
from ~3. 5 MeV to 9 MeV. 

Although copper is not an ideal detector because of its short half-life, it was felt to 
be adequate for purposes of demonstrating the method. Activation fluences were calcu- 
lated for each detector in each flux spectrum. Whenever possible, for actual correla- 
tions using this method, the activation fluences would not be generated in this manner 
but would instead be determined from detector activities and would include the effects of 
irradiation history and half-life. It was necessary to calculate the activation fluences 
for this test of method because this particular set of detectors was not used with the 
existing data. It was felt that calculated activations would be adequate to test and dem- 
onstrate the method. Measured activations would, however, be better because they do 
not contain uncertainties in spectra and cross sections. 

In anticipation of an irradiation temperature effect, temperature was carried as a 
separate entity. Thus, the variables used were the activation fluences and the temper- 
ature at which the specimens were irradiated. 

Four calculations were made to obtain data correlation by four different models: 

(1) A calculation (ref. 11) was performed using neutron fluence above 1 MeV (as- 
suming a Watt fission spectrum) as the only independent variable. Although it is gener- 
ally recognized that the flux above 1 MeV is not a reliable damage predictor, this model 
was selected as representative of the integral flux models. It was felt that the relative 
accuracy of this model compared to other models would be of interest. This calculation 
is referred to as calculation 1. 

(2) Calculation 1 was rerun with temperature added as the second independent vari- 
able. The purpose of this calculation was to determine the effects of temperature when 
used with an arbitrarily selected model. This is called calculation 2. 

(3) A third calculation was run which used activation fluences and temperature as 
independent variables. This calculation is designated calculation 3. 

(4) The fourth calculation which was performed was the same as calculation 3 ex- 
cept that temperature was dropped as an independent variable. This calculation is 
called calculation 4 and was performed to determine whether the temperature effect was 
independent of the selected model, that is, whether the effect would be about the same 
as it is for calculations 1 and 2. 

The damage function method of reference 10 was also used to infer damage values. 
This made it possible to compare results calculated by the activation fluence method 
with those results obtained by a different contemporary technique. 

The data and values of the independent variables used in these calculations are 
shown in tables I and n. The flux correlations (calculations 1 and 2) did not incorporate 
any nonlinear transformations in the fitting procedure as did the activation fluence cor- 
relation (calculation 3). The nonlinear transformation (raising variables to the 1/4 
power) was a numerical expedient that enabled the matrix inversion to be completed. 
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The transformation which makes the X. values lie between 0 and 1. 0 is linear and 
2 1 

has no effect on R and S. When coupled with the 1/4 power transformation, it im- 
proves the matrix conditioning because the transformed values become more homogene- 
ously distributed between 0 and 1. 0. The 1/4 power transformation does affect the R 2 
and S values somewhat. To examine this effect, calculations 1 and 2 were rerun with 
the 1/4 power transformation. The effect on R 2 and S was small; thus, this trans- 
formation did not account for the better agreement of calculation 3, and the conclusions 
remain valid. All results of calculations 1 and 2 shown in this report pertain to those 
calculations in which no 1/4 power transformation was used. 


RESULTS 


In calculations 3 and 4, a high inter correlation was found between the neptunium 
and titanium activities due to the spectral shapes being similar above 1 MeV for about 
70 percent of the data. This resulted in difficulty in inverting the least-squares matrix. 
Since this problem was inherent in the data and because several attempts at transform- 
ing the basic independent variables to avoid excessive covariance proved unsatisfactory, 
the titanium activation was dropped and the neptunium was relied on to sense the com- 
plete fast energy range. 

The least-squares analysis provides for sequential remodeling by dropping less 
significant terms from the equation and refitting until all terms retained are significant 
at some specified level of rejection (ref. 11). The rejection level selected for this study 
was 90 percent. 

The accuracy parameters associated with calculations 1 to 4 are given in table III. 

The parameters include 
o 

(1) R , the square of the multiple correlation coefficient, which is that fraction of 

the total variance in the measured data that is accounted for by the regression 
equation 

(2) S, the standard error of estimate, that is, the standard deviation of fit at the 

mean of the independent variables 

(3) F, the ratio of regression mean sum of squares to residual mean sum of squares 
Comparing the calculated value of F ( F ca i cu i a t e d) to the value of F from the statistical 
tables ( ^tabulated) S ives an indication as to whether the regression equation could have 
the overall accuracy indicated by R 2 by chance alone. A value of F ca i cu i a t e d wel1 
above the tabulated F value indicates only a 5 percent chance that the fit occurred ac- 
cidentally. A high ratio of F calculate d to tabulated is Prerequisite to trusting the 
R 2 , S, and polynomial coefficients. 

(4) F', the ratio of lack-of-fit mean sum of squares to replication mean sum of 

squares 
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When ^calculated is higher than Ftabulated’ the error is probably caused by lack of 
fit and not experimental error. When the reverse is true, the error is probably caused 
by experimental error rather than by the model selection. For a more detailed discus- 
sion of these statistical parameters, see reference 11. 

The fitting coefficients are given for calculation 3 in table IV. Table V is a com- 
parison of measured and predicted points using the first three regression models and 
the selected damage function technique. The following section gives a detailed compari- 
son of the damage function and activation fluence methods. Table VI shows the neutron 
flux distributions that were used to generate activation fluences; and table VII is a sum- 
mary of cross sections obtained for the group structure used by plotting activation 
cross sections (ref. 14) and numerically averaging these over the indicated groups. 

The following are two interpretations that may be made from the data of tables 
m and V: 

(1) These calculations confirm a strong dependence of radiation damage in A302-B 
on the irradiation test temperatures (from 327 to 561 K). This is evident from the large 
increase in the multiple correlation coefficient when the temperature is added to calcu- 
lation 1 (see calculation 2) or calculation 4 (see calculation 3). 

(2) The activation fluence model (for this particular set of data which was assumed 

to be typical and adequate to evaluate relative accuracy) correlates the data much better 

than the integral flux (>1 MeV) method. The activation fluence method explains 93 per- 
2 

cent (see R comparison) of the data variability and has the lowest standard error at 

2 

the fitted points. The lower R for the integral flux (>1 MeV) method was expected. 

The negative coefficient for thermal activation fluence (see table IV) suggests a 
damage annealing effect due to thermal neutrons. The other coefficients also could pos- 
sibly be interpreted by using physical arguments. However, no attempt has been made 
to explain the polynomials from a mechanistic viewpoint since this study was intended 
more to illustrate a technique rather than to produce a usable correlation. The actual 
values of the independent variables might have been different from the values constructed 
using calculated fluxes and cross sections. Also, the data were not taken using an 
orthogonal experiment design; thus, the coefficients may be improved with judicious ex- 
periment planning. It is reasonable to assume, however, that the equations generated 
from this study will have similar accuracy if the independent variables are calculated as 
was done here for other ANDT (in A302-B) values of interest. 

COMPARISON OF ACTIVATION FLUENCE METHOD TO 
DAMAGE FUNCTION METHOD 

It is interesting to note some of the similarities and differences between the activa- 
tion fluence and damage function methods of reference 10. The activation fluence method 
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is similar to the damage function approach of reference 10 in that both are fits of data 
to obtain weighting coefficients. For applications in which activations can be obtained 
experimentally, the activation fluence method has the advantage of not requiring any of 
the complex computer codes which must be used to obtain spectrum information for the 
damage function method. This tends to reduce the propagation and magnification of 
error. The regression analysis with the activation fluence method can be used with the 
minimum number of detectors required for accurate data correlation, whereas the 
damage function model requires detailed spectrum information in many energy groups. 
The regression analysis fit permits treatment of irradiation temperature as an explicit 
independent variable, while the damage function approach of reference 10 assumes that 
the effects of temperature are constant over a temperature range of several hundred 
degrees. A fundamental difference between the regression analysis method and the 
damage function method of reference 10 is that the damage function method requires an 
additional separate correlation between damage and total neutron fluence in order to 
predict the damage that will occur due to an arbitrary spectrum and fluence. Regres- 
sion analysis using a single correlation predicts damage within the limits of the data 
used to derive the coefficients. Finally, the activation fluence method allows explicitly 
for interaction effects between detectors and between detectors and temperature. This 
effect is implicit in the damage function method. 

The activation fluence polynomial can be solved for that total fluence required to 
produce a specified change in nil-ductility temperature by replacing the time-integrated 
activations by the product of total fluence and fluence-weighted activation cross sections, 
provided that these cross sections are available or can be generated. The damage func- 
tion method yields the required total fluence directly. 

The activation fluence method correlates these data as well as the damage function 
method (11. 5 K standard error for the activation fluence method compared to 12. 6 K 
standard error for the damage function method). Table V shows a detailed comparison 
of measured and calculated values. 

CONCLUSIONS AND RECOMMENDATIONS 

Based on the results of this study, it is concluded that the activation fluence method 
may be an accurate and relatively simple method for correlating radiation damage. Its 
disadvantage is that for data having high covariance, the aliasing of effects will prevent 
physical interpretation of individual coefficients. However, this does not hamper its 
predictive capability, which is the goal of a regression analysis of unplanned data. 

It is also concluded that regression analysis provides an easy way to accommodate 
temperature as an additional variable. Finally, it is noteworthy that a three-energy- 
group spectrum description appears adequate for correlating A302-B nil-ductility data. 
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This work in applying the activation fluence method points out the need for standard- 
ization of the method of determining the spectrum variable. While the number and 
choice of detectors will depend on the property and temperature range for the material 
being studied, a complete description of the dosimetry method will permit other inves- 
tigators to use the data generated in exploring other methods of data correlation. The 
methods used at the Plum Brook Reactor Facility are described in part V of "Standard 
Guides to the Design of Experiments for the Plum Brook Reactor Facility. " 

With the activation fluence method, when detectors can be placed in-pile with the 
specimens, they may be left in-pile for the duration of the irradiation if long half-life 
detectors are selected exclusively. If some short half-life detectors are used, at least 
one long half-life detector should be used so that the detectors can be normalized to 
comply with the definition of activation fluence given in equation (2). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 24, 1971, 

112-27. 
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TABLE I. - SUMMARY OF DATA ANALYZED 





H 

r T /no MeV 

4 A Mev ®( E ' dt 

(b) 

/•T /'10 MeV 

^ J Q <p(E, t)o(E)dE dt 

(cobalt) 

/• T r 10 MeV 

jg y Q <p(E, t)o(E)dE dt 

(copper) 

/■T /-lOMeV 

J Q J q <p{E, tHE)dE dt 

(neptunium) 

Irradia- 
tion tem- 
perature, 
K 

1 

LITR 

C- 

18 

77. 8 

0. 50X10 19 

21. 4X10 19 

0. 207X10 19 

0. 950x10 19 

505 

1 



C- 

18 

77.8 

. 50 

21.4 

.207 

.950 

506 

2 



C- 

35 

86. 1 

3. 1 

164.5 

1.543 

6.694 

561 

2 



C- 

55 

94.4 

3. 1 

164.5 

1. 543 

6.694 

561 

3 



C- 

18 

94. 4 

. 50 

21. 4 

.207 

. 950 

400 

4 





72.2 

. 50 

21.4 

.207 

.950 

477 

5 





111. 1 

.60 

25.7 

. 249 

1. 140 

366 

6 





141.7 

1.8 

77.0 

. 747 

3. 420 

327 

7 



C- 

28 

122.2 

1.24 

72,0 

.356 

2. 279 



8 



C- 

53 

119.4 

1.2 

74. 2 

. 586 

l. 422 



9 





163.9 

1. 94 

119. 9 

.948 

2. 299 



10 





175.0 

2.00 

123.6 

. 977 

2. 370 



11 





183. 3 

3. 50 

216.4 

1. 710 

4. 148 



12 



C- 

55 

127.8 

.85 

451.2 

. 423 

1. 835 



13 



C- 

49 

113.9 

.75 

82.2 

.663 

1. 801 



14 

IRL 

11. 75 cm (4. 625 in. ) 

58.3 

.26 

4. 1 

.048 

. 339 



15 



14. 29 cm 

(5. 625 in. ) 

44. 4 

. 16 

1.6 

. 041 

.274 



16 



16. 83 cm 

(6.625 in.) 

27.8 

. 10 

1.0 

. 035 

.211 



17 



19. 38 cm (7. 625 in. ) 

27.8 

.06 

. 8 

! .029 

. 158 



18 



21. 91 cm (8. 625 in. ) 

19.4 

.04 

.8 

r . 027 

. 127 



19 

BGR 

W 

-44 

113.9 

. 55 

4. 1 

. 130 

1. 572 



20 

BGR 

W 

-44 

133.3 

.75 

5. 5 

. 177 

2. 144 



21 

HWCTR 

Gray rod 

105.6 

.73 

1030. 0 

1. 199 

1. 515 



22 

LITR 

C 

55 

111. 1 

1. 4 

74. 3 

.698 

3.023 

527 

23 

Yankee 

Vessel wall 

61. 1 

. 22 

19. 1 

.041 

. 221 

544 

24 



Accelerated 

125.0 

5.0 

1855.6 

. 965 

6. 306 

547 

25 



Accelerated 

144. 4 

7.0 

2597. 5 

1.351 

8. 831 

547 

26 



Accelerated 

I 172.2 

9.0 

3339. 4 

1.736 

11. 350 

541 

27 

LITR 

C 

-55 

86. 1 

1.5 

79.6 

. 747 

3. 238 



28 



C 

-18 

86. 1 

3.3 

141. 1 

1. 369 

6.268 



29 



C 

-18 

91.7 

3.0 

128.3 

1. 244 

5.700 



30 



C 

-55 

105.6 

4.7 

249. 5 

2.341 

J 

10. 144 




a See refs. 12 and 13 for a description of the exact position of the specimen relative to the reactor core. 

Watt fission spectrum is assumed for this integral. All other integrals in this table use the calculated (assumed actual) flux distribution. 






TABLE n. - MODEL DESCRIPTION 


Calcula- 

tion 


Primary model 


Transformation 


a Y = A fl + AjZ + A n Z 2 


r IQ /” r Z’ 10 MeV 1 

Z = [lO y iMeV <p(E> t)dE dtj 


Y - a q + A 1 G 1 + A 2 G 2 + A 12 G 1 G 2 
+ A 11 G 1 + A 22 G 2 


/~T /"lO Mev 7 

10 ~ 19 4 -4m eV <?>(E,t)dEdtJ 

G 2 = Irradiation temperature 


Y - A 0 + AjXj + A 2 X 2 + A 3 X 3 


+ A 4 X 4 + A 12 X 1 X 2 + A 13 X 1 X 3 


+ A 14 X 1 X 4 + A 23 X 2 X 3 


+ A 24 X 2 X 4 + A 34 X 3 X 4 + A 44 X 4 


■10 MeV 1 


ip(E, t)o(E)dE dt 

Co- 59 

•10 MeV t 


ip(E, t)o(E)dE dt 

Cu-63 

10 MeV "| 


<p( E, t)u(E)dE dt 

Np-237 


Y - A 0 + AjXj + A 2 X 2 + AjXg 


+ a 12 x 1 X 2 + Aj^XjXj 


+ A 23 X 2 X 3 


a . - Irradiation temperature 

1/41 


X i = 


/ 

min \ 

/ O.J- 


l max 

1 

mini 


a Y is the change in nil-ductility transition temperature for A302-B ferritic steel, i. e. , ANDT of table I. 

b^max and o min are the largest and smallest activation fluences (for the i th detector) that were generated 
11 4 

from the data used in the regression analysis. The linear transformation from a. to X ; codes the data 

so that it will be between 0 and 1. 0. 


TABLE III. - SUMMARY OF ACCURACY PARAMETERS 


Statistical parameter 

(a) 

Symbol 

Calculation 1^ 

Calculation 2 

Calculation 3 

Calculation 4 


Primary 

model 

Reduced 

model 0 

Primary 

model 

Reduced 

model 

Primary 

model 

Reduced 

model 

Primary 

model 

Reduced 

model 

Regression sum of squares 
Total sum of squares 

R 2 

0. 355 

0. 302 

0. 788 

0. 764 

0. 931 

0. 925 

0.644 

0. 604 

(Residual mean square)^/ 2 

S 

35. 0 

35.8 

21. 2 

21. 2 

13. 8 

13. 1 

28. 0 

27. 9 

Regression mean square 

F 1 

.calculated 

7. 97 

13.6 

19. 3 

46. 8 

24. 6 

38. 4 

7. 53 

14. 1 

Residual mean square 

d F 

tabulated 

3. 33 

4. 18 

2. 59 

3. 37 

2. 32 

2. 51 

2. 49 

2. 99 

Lack-of-fit mean square 

p' 

calculated 

75. 9 

79. 2 

28. 0 

27. 7 

12. 1 

10. 7 

49. 1 

48. 3 

Replication mean square 

d F' 

tabulated 

19. 5 

19. 5 

19. 5 

19. 5 

19. 4 

19. 5 

19. 5 

19. 5 


a Total sum of squares is equal to regression sum of squares plus residual sum of squares. Residual sum of squares is 


equal to lack-of-fit sum of squares plus replication sum of squares. 

^See RESULTS section of text for a description of each calculation. 

c Reduced model refers to model obtained after all terms not significant at 90 percent level have been deleted by stepwise 
backward elimination procedure. 

^Probability value is 95 percent. 
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TABLE IV. - SUMMARY OF ACTIVATION 


FLUENCE MODEL COEFFICIENTS 
FOR CALCULATION 3 

[aNDT (in K, primary model) = A Q + AjXj + A 2 X 2 

AgXg + A 4 X 4 + A^X^ + AjgXjXg 

+ a 14 X i X 4 + A 23 X 2 X 3 + A 24 X 2 X 4 + A 34 X 3 X 4 
+ A 44 X 4 ; ANDT (in K, reduced model) = Aq 
+ A 1 X 1 + AgXg + A 4 X 4 + A 14 XjX 4 + A 2 gX 2 X 3 

+ A 34 X 3 X 4 + A 44 X 4-] 


Coefficient or 
constant 

Primary model 

Reduced model 

A 0 

-8. 33 

-8. 09 

A 1 

-177.8 

-73. 3 

A 2 

105. 4 


A 3 

87.6 

129.3 

A 4 

241.5 

265. 6 

A 12 

-120.0 


A 13 

241. 3 


A 14 

243. 9 

190. 2 

A 23 

88. 1 

180.0 

A 24 

461. 1 


A 34 

-850.0 

-343. 4 

A 44 

-139.0 

-169. 3 


a See table n for definitions of the independent var- 


iables X. 




TABLE V. - COMPARISON OF PREDICTED AND MEASURED DATA 


Measured 
change in 
nil-ductility 
temperature, 
ANDT, 

K 

(see table I) 

Predicted change in nil-ductility temperature ANDT, K; 

Measured - Calculated ^qq 
M easured 

and error, 


Calculation 1 

Calculation 2 

Calculation 3 

Damage function 

Calculated 

Percent 

error 

Calculated 

Percent 

error 

Calculated 

Percent 

error 

Calculated 

Percent 

error 

77.8 

87. 3 

-12. 3 

71. 9 

+7.6 

71. 4 

+8.2 

88. 

3 

-13.6 

77. 8 

87. 3 

-12. 3 

71. 9 

+7.6 

71. 4 

+8. 2 

88 

3 

-13.6 

86. 1 

115. 8 

-34. 5 

93. 1 

-8. 1 

97. 7 

-13. 5 

(a) 


94. 4 

115.8 

-22.6 

93. 1 

+4. 1 

97.7 

-3. 5 

(a) 


94. 4 

87. 3 

+7. 5 

80.2 

+ 15. 1 

98.6 

-4. 4 

88 

3 

+6. 5 

72. 2 

87. 3 

-20. 1 

74. 1 

-2. 6 

78. 1 

-8.2 

88 

3 

-22. 3 

111. 1 

88. 4 

+20.4 

86.8 

+21.9 

112.2 

-1.0 

95 

6 

+ 14.0 

141.7 

101.6 

+28.3 

146.2 

-3. 2 

161.0 

-13. 6 

138 

3 

+2. 4 

122. 2 

95. 4 

+21. 9 

120. 2 

+1.7 

127. 1 

-4.0 

122 

2 

+0.0 

119. 4 

95.0 

+20. 5 

118. 3 

+.9 

118. 3 

+.9 

125 

6 

-5. 1 

163. 9 

103. 1 

+37. 1 

152. 7 

+6.8 

145.6 

+ 11. 2 

144 

4 

+ 11. 9 

175. 0 

103. 8 

+40. 7 

155. 4 

+11. 2 

147. 4 

+15. 8 

145 

6 

+ 16. 8 

183. 3 

120. 2 

+34. 4 

225. 1 

-22. 8 

187. 4 

-2. 2 

167 

8 

+8. 5 

127. 8 

91. 2 

+28. 7 

102. 1 

+20. 1 

124. 8 

+2. 3 

113 

9 

+ 10. 9 

113. 9 

90. 1 

+20.9 

97. 4 

+ 14. 4 

129. 5 

-13. 7 

114 

4 

-. 5 

58. 3 

84. 7 

-4J>. 2 

74. 7 

-28. 1 

54. 8 

+6.0 

46 

1 

+21. 0 

44. 4 

83. 6 

-88. 1 

70. 1 

-57.6 

50. 7 

-14. 1 

37 

8 

+ 15.0 

27. 8 

82. 9 

-199.6 

67. 3 

-142.2 

43. 8 

-57.6 

29 

4 

-6. 0 

27. 8 

82. 5 

-197. 0 

65.4 

-135. 5 

36. 4 

-31. 9 

16 

7 

+40. 0 

19. 4 

82. 3 

-323. 2 

64. 5 

-231. 7 

4.0 

+79. 4 

9 

4 

+51. 4 

113. 9 

87. 9 

+22. 8 

88. 2 

+22. 6 

109. 7 

+3. 7 

106 

7 

+6. 3 

133.3 

90. 1 

+32. 4 

97.4 

+26. 9 

124. 1 

+6.9 

118 

9 

+ 10.8 

105. 6 

89. 8 

+14. 9 

96. 5 

+8. 6 

112. 5 

-6. 5 

107 

.2 

-1.6 

111. 1 

97. 2 

+12. 5 

83.7 

+24. 7 

83. 8 

+24. 6 

(a) 


61. 1 

84. 3 

-37. 9 

65. 4 

-7.0 

72.7 

-19.0 




125.0 

136. 7 

-9. 3 

122.7 

+1. 9 

131. 6 

-5. 3 




144. 4 

158. 6 

-9.8 

146.7 

-1. 5 

147. 3 

-2. 0 




172. 2 

180. 5 

-4. 8 

170. 7 

+. 9 

161.3 

+6. 3 




86. 1 

98. 3 

-14. 1 

77. 4 

+10. 2 

76. 3 

+ 11. 4 




86. 1 

118.0 

-37. 0 

95. 1 

-10. 4 

92. 2 

-7. 1 




91. 7 

114. 7 

-25. 2 

92. 1 

-. 5 

88. 9 

+3. 0 




105.6 

133. 4 

-26.3 

108. 8 

-3. 1 

116.2 

-10. 0 




Standard 
error, K 

35.8 

21. 2 

13. 1 
b 11.5 

12.6 


a Cannot be calculated using the damage function because irradiation temperature is greater than 
506 K (450° F). 

b For same sample of 21 points used for damage function. 
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C-18 C-28 C-49 


7. 79 MeV 
6. 07 MeV 
4. 72 MeV 
3. 68 MeV 
2. 87 MeV 
2. 23 MeV 
1. 74 MeV 
1. 35 MeV 
1.05 MeV 
. 821 MeV 
. 639 MeV 
. 498 MeV 
. 388 MeV 
. 302 MeV 
. 235 MeV 
. 183 MeV 
. 414 eV 


TABLE VI. - SUMMARY OF RELATIVE NEUTRON FLUX USED 


Reactor 


IRL 


Test specimen location 


Relative neutron flux 


BGR HWCTR 


11.75 cm 

14. 29 cm 

16. 83 cm 

19. 38 cm 

21. 91 cm 

W-44 

(4.625 in.) 

(5.625 in.) 

(6. 625 in. ) 

(7.625 in.) 

(8. 625 in. ) 



Gray rod Pres- 
sure 


0. 00257 
. 0966 
.267 

. 401 
. 565 
1.043 

1. 444 
1.757 
L. 834 
1. 863 
1.944 
1. 902 
1. 781 
1. 766 
1.686 
1.638 
9.646 
8 


The group fluxes are the fluxes integrated over each group. All nonthermal flux distributions are calculated distributions at the specimen locations, 
measured at these locations. 


fluxes were 



































TABLE VII. - SUMMARY OF SELECTED CROSS SECTIONS 


Energy 

group 

Lowest 

energy 

Cross sections, barns 

Co 59 (n,y)Co 60 

a Cu 63 (n, y)Cu 64 

Np 237 (n, f)N fp 

i 

7. 79 MeV 

0. 00012 

0. 0034 

2. 35 

2 

6. 07 MeV 

. 00014 

. 0037 

1. 80 

3 

4. 72 MeV 

. 00021 

.0042 

1. 42 

4 

3. 68 MeV 

. 00045 

. 0048 

1. 50 

5 

2.87 MeV 

. 00095 

.0056 

1. 55 

6 

2. 23 MeV 

.0016 

. 0066 

1.62 

7 

1. 74 MeV 

.0025 

. 0075 

1. 67 

8 

1. 35 MeV 

. 0031 

.0083 

1.64 

9 

1. 05 MeV 

.0039 

. 0094 

1. 58 

10 

. 821 MeV 

.0044 

.0107 

1.33 

11 

. 639 MeV 

.0042 

.0125 

1. 00 

12 

. 498 MeV 

.0040 

.0145 

.650 

13 

. 388 MeV 

.0053 

. 0170 

.270 

14 

. 302 MeV 

.0062 

.0210 

. 110 

15 

. 235 MeV 

.0065 

.0235 

.0400 

16 

. 183 MeV 

. 0068 

.0280 

.0275 

17 

. 414 eV 

3.61 

. 153 

.0109 

18 

0 

37 

0 

.019 


a These cross sections are with the standard 0. 0508 cm cadmium 
cover. 
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